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Stress  generation  due  to  Li  ion  insertion  into/extraction  from  LiMn204  particles  is  studied  with  a  mathe¬ 
matical  model  for  a  lithium  ion  battery  with  pure  LiMn204  or  mixed  LiMn204  and  LiNi0.sCo0.15Al0.05O2 
cathode.  The  simulated  stress  profile  in  a  pure  LiMn204  electrode  shows  nonuniformity  across  the  positive 
electrode.  The  cathode  blended  model  predicts  that  the  stress  generated  in  the  LiMn204  particles  is 
reduced  at  the  end  of  discharge  due  to  adding  LiNio.sCoo.15Alo.05O2  to  the  cathode.  The  effect  of  the  variation 
in  the  blend  ratio  on  the  stress  generation  is  also  investigated. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Spinel  LiMn204  (LMO)  is  one  of  the  most  attractive  positive 
electrodes  for  high  power  applications  such  as  hybrid/plug-in 
hybrid  electrical  vehicles  (HEV/PHEVs).  It  shows  great  benefits  in 
the  area  of  cost,  safety  and  power  density,  which  are  major  con¬ 
cerns  of  large  format  applications  in  the  automotive  industry. 
However,  this  spinel  material  exhibits  serious  capacity  fade  during 
cycling  or  storage  at  elevated  temperature  which  makes  it  less 
competitive  with  other  cathode  materials  [1].  Recent  studies  have 
revealed  that  capacity  fade  in  spinel  LMO  cathodes  is  mostly  due  to 
manganese  (Mn)  dissolution  and  its  toxic  effect  on  the  carbon 
anode  [2-4].  There  are  at  least  two  possible  reasons  for  the  Mn 
dissolution:  (1)  acid  attack  and  a  disproportionation  reaction  of 
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Mn(III)  on  the  particle  surface  [5-7];  (2)  the  instability  of  the  two- 
phase  structure  in  the  charged  state  which  leads  to  the  loss  of  MnO 
and  the  dissolution  of  Mn(II)  to  form  a  more  stable  single-phase 
structure  [8]. 

Several  improvements  have  been  made  to  increase  the  long 
term  stability  of  the  LMO  based  cells  by  using  surface  treatment  or 
element  substitution.  For  example,  Amatucci  et  al.  [9]  have  inves¬ 
tigated  the  effects  of  modifying  the  surface  of  LMO  on  the  cell 
performance  at  high  temperatures.  The  surface  improvement  helps 
to  reduce  Mn  dissolution  and  retain  capacity.  Doped  spinel 
LiMn2_xMx04  (where  M  =  Li,  Co,  Cr,  Ni  etc.)  materials  have  also 
been  reported  to  be  effective  for  the  improvement  of  the  cell  high 
temperature  performance  [10-13].  The  benefit  of  doping  elements 
in  LMO  has  been  attributed  to  the  stabilization  of  the  spinel 
structure.  Among  them,  the  doped  material  with  lithium  substi¬ 
tution  of  Mn  sites  has  higher  theoretical  capacity  than  others 
because  of  the  light  element  lithium,  but  the  doped  material  retain 
the  native  low  capacity  of  LMO.  In  order  to  increase  the  reversible 


366 


Y.  Dai  et  al.  /  Journal  of  Power  Sources  247  (2014)  365-376 


capacity  after  the  cation  doping  while  keeping  the  other  electro¬ 
chemical  performance,  F-  ions  were  introduced  in  place  of  O2-  ions 
[14,15].  Unfortunately,  the  cost  for  the  production  of  the  doped 
materials  is  significantly  high. 

A  promising  approach  was  first  provided  by  Numata  et  al.  to 
improve  the  capacity  retention  for  LMO  based  electrodes  by  adding 
LiNio.8Coo.2O2  [16].  They  found  that  both  the  Mn  dissolution  from 
the  LMO  and  the  acid  generation  decreased  with  the  addition  of 
LiNio.8Coo.2O2.  Myung  et  al.  [17  investigated  the  behavior  of  the 
electrode  with  the  LMO  and  LiNio.8Coo.15Alo.05O2  (NCA)  mixture. 
They  found  that  a  cell  with  the  mixed  electrode  had  better  cycle 
performance  at  elevated  temperature  than  a  cell  with  a  pure  LMO 
electrode.  In  addition,  the  capacity  of  the  LMO  based  electrode  has 
been  enhanced  by  adding  high  capacity  NCA.  Recently,  Tran  et  al. 

[18]  reported  on  the  electrochemical  and  thermal  behavior  of  LMO 
and  NCA  blended  electrodes.  Their  results  indicated  that  the  NCA 
addition  increased  the  capacity  of  blends  and  reduced  the  Mn 
dissolution  from  the  spinel,  but  harmed  the  rate  capacity;  On  the 
other  hand,  with  the  existence  of  LMO  the  heat  generation  from  the 
blends  is  less  than  that  from  the  pure  NCA.  According  to  their  study, 
blends  with  33.3  wt%  NCA  have  the  best  behavior.  Manthiram  et  al. 

[19]  studied  the  suppression  of  the  Mn  dissolution  in  spinel  cath¬ 
odes  by  mixing  LMO  or  its  doping  derivatives  with  layered  oxide 
cathodes  such  as:  LiCo02  and  LiNio.85Coo.15O2.  They  observed  that 
LiCo02  was  more  effective  than  LiNio.85Coo.15O2  in  improving  the 
storage  and  cycle  performance  of  LMO  based  electrodes.  Smith  et  al. 

[20]  showed  that  the  addition  of  Li[Nii/3Mni/3Coi/3]02  to  the  LMO 
electrode  produced  a  great  improvement  in  the  capacity  retention. 
Also,  it  was  found  that  the  Li[Nii/3Mni/3Coi/3]02  helped  to  reduce 
the  Mn  dissolution  from  the  LMO  electrode.  Now,  Li-ion  batteries 
with  LMO  mixed  with  other  insertion  active  material  as  cathodes 
have  been  commercialized  such  as  Samsung  SDI  Li-ion  batteries. 

The  mechanical  degradation  under  cycling  is  one  of  the  most 
important  failure  mechanisms  of  modern  Li-ion  batteries.  It  has 
been  reported  that  the  fractures  and  cracks  exist  in  several  elec¬ 
trode  materials  after  Li  insertion/extraction  [21-24].  For  example, 
Lim  et  al.  [23]  have  observed  an  evident  disruption  in  the  LMO 
particle  after  electrochemical  cycling.  There  is  about  6.5%  volume 
change  when  the  LMO  is  lithiated  from  empty  state  to  full  state 
[25].  Flowever,  the  particle  volume  change  inside  particle  is  not 
uniform.  This  gives  rise  to  stress  and  produces  cracks  or  fractures  in 
the  particle.  Furthermore,  the  particle  fracture  will  generate  fresh 
active  material  surface  on  which  the  new  solid  electrolyte  interface 
(SEI)  film  is  formed.  In  the  meanwhile,  the  particle  fracture  can  also 
lead  to  loss  of  connection  of  the  active  material  particle  to  the 
conductor  [26]. 

Mathematical  models  have  been  largely  applied  to  investigate 
the  stress  generation  and  the  effect  of  stress  on  battery  degrada¬ 
tion.  Christensen  et  al.  [27,28]  developed  an  integrated  model  based 
on  the  theories  of  transport  and  elasticity  to  simulate  the  stress 
inside  active  material  particles.  Zhang  et  al.  [29]  studied  the 
intercalation-induced  stress  in  LMO  particles  by  considering  the 
intercalation-induced  stress  as  being  analogous  to  thermal  stress. 
Their  approach  showed  similar  results  to  those  of  Christensen  et  al. 
Without  considering  the  stress-driven  diffusion  in  Zhang  et  al.’s 
approach,  Cheng  et  al.  [30]  obtained  an  analytical  solution  for  the 
generation  of  stress  in  a  spherical  particle  under  galvanostatic  and 
potentiostatic  operation. 

Unlike  the  studies  above  on  single  particles  only,  there  are 
several  papers  on  the  stress  generation  in  a  porous  composite 
electrode.  Garcia  et  al.  [31]  developed  a  two-dimensional  model 
with  particle  distribution  to  compute  stress  in  the  porous  electrode. 
Also,  their  model  account  for  the  potential  and  Li-ion  concentration 
distribution  in  the  solution.  However,  their  model  is  very  compu¬ 
tationally  intensive  and  the  stress-driven  diffusion  has  been 


ignored.  Golmon  et  al.  [32  combined  the  mechanical  equation  and 
porous  electrode  theory  (the  pseudo  two  dimensional  (P2D) 
model)  together  to  investigate  the  stress  distribution  in  porous 
electrode.  They  applied  homogenization  techniques  to  relate  the 
parameters  between  the  micro-scale  and  the  macro-scale.  Simi¬ 
larly,  Renganathan  et  al.  [33]  studied  the  stress  generation  in  a 
LiCo02/Carbon  battery  including  the  phase  change  in  the  LiCo02. 
Christensen  [34  simulated  the  stress  generation  in  a  porous  elec¬ 
trode  by  incorporated  the  model  for  stress  generation  in  single 
particle  [27,28]  to  the  P2D  model.  Recently,  Bower  et  al.  [35  pro¬ 
posed  a  finite-strain  elastic-plastic  model  to  predict  the  variation  of 
stress  and  electric  potential  in  a  one-dimensional  half-cell. 

However,  the  stress  generation  in  an  electrode  with  multiple 
active  materials  has  not  been  published  in  the  open  literature.  To 
obtain  a  better  understanding  of  the  properties  of  Li-ion  batteries 
with  blended  material  cathodes,  we  present  a  theoretical  study  of 
the  mechanical  behavior  of  electrodes  made  from  mixing  LMO  and 
NCA.  The  mechanical  equations  capturing  the  stress  generation  in 
the  spherical  particles  in  the  electrodes  are  incorporated  into  the 
classical  P2D  model.  To  elucidate  the  effect  of  mixing,  simulations 
are  conducted  and  compared  for  both  pure  LMO  electrodes  and 
blended  electrodes  with  LMO  and  NCA.  Simulation  results  for  a 
mixed  electrode  show  that  the  stress  generated  in  the  LMO  parti¬ 
cles  is  significantly  reduced  at  the  end  of  discharge  due  to  adding 
NCA  to  the  electrode.  This  stress  reduction  might  partially  explain 
why  experimental  investigations  show  that  the  performance  of  the 
LMO  based  electrodes  can  be  improved  by  adding  NCA.  We  also 
investigated  the  mixing  ratio  on  the  stress  generation  of  the 
blended  electrode.  Simulations  for  mixed  electrodes  with  different 
blend  ratios  show  that  the  reduced  stress  region  increases  by 
increasing  the  NCA  ratio. 

2.  Model  development 

In  this  study,  we  have  combined  P2D  model  and  mechanical 
equations  to  study  the  stress  generated  due  to  Li  insertion  into / 
extraction  from  the  active  material  particles  in  a  lithium  ion  bat¬ 
tery.  This  approach  is  similar  to  Renganathan  et  al.’s  work  to  study 
the  LiCo02/Carbon  battery  [33],  but  we  have  extended  their  work  to 
electrodes  with  multiple  active  materials  in  the  cathode.  Further¬ 
more,  because  no  phase  change  is  considered  during  Li  insertion 
into/extraction  from  the  active  material  particles,  we  could  express 
the  stress  profile  as  function  of  concentration.  To  obtain  the  stress 
profile  in  the  P2D  model,  we  modified  the  diffusion  equation  in  the 
classic  P2D  model,  which  is  shown  in  detail  below. 

The  Li-ion  battery  considered  in  this  study  consists  of  three 
regions:  a  positive  electrode  (pure  LMO  or  LMO  mixed  with  NCA),  a 
separator,  and  a  carbon  negative  electrode  as  shown  in  Fig.  1.  For  an 
electrode  that  has  more  than  one  active  material,  it  is  assumed  that 
the  materials  in  the  electrode  are  well  mixed  and  their  properties 
(such  as:  diffusion  coefficients,  reaction  rate  constants  and  so  on) 
remain  the  same  as  those  in  a  pure  material  electrode.  The  pa¬ 
rameters  for  the  mixed  positive  electrode  used  in  this  model  were 
obtained  from  the  work  by  Albertus  et  al.  [36],  in  which  the  elec¬ 
trochemical  behavior  of  a  mixed  cathode  half-cell  was  simulated.  In 
the  following,  we  will  briefly  introduce  the  mechanical  equations 
and  the  porous  electrode  model  equations. 

2.1.  Mechanical  stress  model 

In  this  work,  the  stress  generated  due  to  Li-ion  intercalation / 
deintercalation  process  is  calculated  using  an  approach  similar  to 
that  due  to  volumetric  thermal  expansion.  The  thermal  stress 
analysis  of  an  isotropic  media  can  be  described  by  Hooke’s  law  with 
additional  thermal  expansion  term  [37]: 
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Fig.  1.  Schematic  of  a  Li-ion  battery. 
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where  e- y  and  ay  represent  the  strain  and  stress  tensor  components, 
respectively;  E  is  Young’s  modulus,  v  is  Poisson’s  ratio,  5y  is  Kro- 
necker  delta,  a  is  the  thermal  expansion  coefficient,  and  AT  repre¬ 
sents  the  temperature  difference  from  the  original  value. 

Analogously,  in  the  case  of  small  deformation,  a  stress-strain 
relationship  including  a  concentration  diffusion  in  an  isotropic 
media  can  be  represented  as  38]: 

£ij  =  £  [(1  +  v)aij  ~  v(JklAj]  3  ^  ^ij  (2) 

where  Q  is  the  partial  molar  volume,  Acs  =  cs  -  c°s  is  the  concen¬ 
tration  change  of  the  diffusing  species  from  original  value,  c°s  and  cs 
is  the  concentration  of  the  diffusing  species. 

For  a  spherical  particle,  there  are  only  two  independent  com¬ 
ponents  in  the  stress  tensor:  the  radial  component  (o>)  and  the 
tangential  component  (oy).  Therefore,  the  equilibrium  condition  for 
the  stress  tensor  in  a  particle  is  simplified  as  follows  [29]: 


^L  +  l{ar-at)  =  0  (3) 

Also,  strain-stress  relations  can  be  written  according  to  equa¬ 
tion  for  a  spherical  particle  [29], 

£r  =  -£  (0>  —  2v(Tt)  +  jAcs 

et  =  l  fit  ~  v{at  +  <7r)]  +  § Acs 

The  stresses  can  be  written  as  a  function  of  strains 

ar  =  (i+„)fl-2„)  [C  -  v>r  +  2 vet  “  (1  +  v) 

at  =  (l+i/)fl-2i/)  [V£r  +  et  -  (1  +  v) 

Also,  strains  can  be  related  to  the  displacement  u  in  a  spherical 
particle  as  follows: 


P  —  du 
£r  ~  df 


£t  —  7 


(6) 


The  u  is  defined  as  the  difference  between  the  radial  position  of 
a  lattice  at  time  t,  and  the  initial  radial  position  of  that  lattice 
element. 

By  substituting  Eq.  (6)  into  (5),  the  stress  tensor  can  be 
expressed  in  terms  of  displacement  u  as  follows: 

T^r)[0-^  +  ^UT-0+v)¥\ 

Substitute  Eq.  (7)  into  the  equilibrium  condition  (Eq.  (3))  yields 
the  displacement  equation  which  is  shown  as  following: 


d2u  2du  2u_l+y£dcs 
dr2  +  r  dr  r2  -  1  -  v  3  dr 

Integrating  Eq.  (8)  yields: 


du  _u  1  +vQ  _ 

-j-  +  2-  =  ~z - -C-\~Z\ 

dr  r  1  -  v  3 

Integrating  Eq.  (9)  yields: 

1  (  [\+vQ  2a  zir3  7 
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(9) 


(10) 


where  Z\  and  Z2  are  integration  constants.  These  two  constants  can 
be  obtained  from  the  two  boundary  conditions:  the  radial  stress  is 
zero  on  the  particle  surface  (r  =  R0) 

&r\(r  =  R0)  =  0  (11) 

and  the  stress  remains  finite  at  the  particle  center  (r  =  0),  that  is: 

=  0  (12) 


dar 

dr 
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Substituting  the  solution  for  u  into  Eq.  (7),  results  in  the 
following  expressions  for  the  two  stress  components  [29]: 
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where  £  is  a  dummy  integration  variable. 

Finally,  the  hydrostatic  stress,  ah,  is  defined  and  given  as  follows 

[29,38]: 


°h 


(Jr  +  2  at 
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2  QE 

90^7) 
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csr2dr  -  cs 


(14) 


Because  the  Eqs.  (13)  and  (14)  are  derived  without  assuming 
that  there  is  constant  applied  current/flux  on  the  particle  surface,  it 
is  reasonable  to  combined  this  with  the  P2D  model  as  shown  below. 


2.2.  Porous  electrode  model 

The  porous  electrode  model  (P2D  model)  has  been  applied  to 
simulate  the  behavior  of  Li-ion  batteries  in  many  literature  39-41  ]. 
The  mass  balance  and  the  charge  balance  are  considered  in  the  both 
solution  phase  and  the  solid  phase.  We  only  briefly  introduce  the 
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model  here  as  follows.  A  detailed  explanation  of  the  porous  elec¬ 
trode  models  can  be  found  elsewhere  [39-41  ]. 

The  mass  balance  for  Li+  in  the  electrolyte  is  given  as  follows: 

e2j^f  =  -v(-DeffLi+J.VcLi+)  +  1  Ft+/j,  j  =  P,s,n  (15) 

where  £2j  is  the  porosity  in  cell  component  j  (j  =  p  for  the  cathode, 
j  =  s  for  the  separator,  and  j  =  n  for  the  anode),  cLi+  is  the  con¬ 
centration  of  Li+,  Deff  Li+J-  is  the  effective  diffusivity  of  Li+  in  the 
electrolyte,  t+  is  the  transport  number  of  Li+,  Fis  Faraday’s  constant, 
and  Ij  is  the  local  total  current  density  (A  m  3).  The  potential  in  the 
solid  phase,  0i,  is  given  by  Ohm’s  law  as  follows: 

V-(  — ffrffjVtfi  )=-ljJ  =  P,n  (16) 


the  surface  of  particle  type  k.  The  effective  properties  are  discussed 
and  presented  in  the  Appendix. 

2.3.  Diffusion  in  the  solid  phase 

The  diffusion  equation  is  modified  by  including  the  stress-driven 
diffusion.  The  pore  wall  flux,  Jk,  for  particles  type  k  including  the 
stress-driven  diffusion  can  be  expressed  as  following  29,38]: 

Jk  =  ~Ds,k  (vcs,k  -  ^V<r Klfj  (23) 

where  Ds  /<  is  the  diffusion  coefficient  in  particles  type  k,  cs,k  is  Li  ion 
concentration  in  particles  type  k,  and  dhj<  is  the  hydrostatic  stress  of 
particle  k.  The  mass  balance  for  Li  ions  in  particles  type  k  is  given 
by: 


where  creffj  is  the  effective  conductivity  in  the  solid  phase. 
The  solution  phase  potential,  is  given  by: 


®cs,k 

~~9 T 


-v-A 


(24) 


(l  +jg:)vlncu,)  (17) 

=  IjJ  =  p,s,n 

where  /ceffj  is  the  effective  conductivity  in  the  solution  phase,  R  is 
the  gas  constant,  T  is  the  temperature,  /  is  the  ionic  activity 
coefficient. 

The  local  total  current  density,  Ij,  is  the  sum  of  current  densities 
for  all  types  of  particles 

h  =  j  =  P-n  (18) 

k 

where  kj  is  particle  types  in  region  j.  In  the  positive  electrode,  there 
are  two  types  of  particles:  LMO  and  NCA;  in  the  negative  electrode, 
there  is  only  one  particle  type,  is  the  specific  surface  area  of 
types  particle  k  and  is  defined  as 


with  two  boundary  conditions  as  following: 

Jk  I  (r  =  0)  =  0 

//<l(r  =  R0,t)  =  T 

Substitute  Eqs.  (23)  into  (24), 

®csj< 


at 


=  -V-  -D, 


Js,k 


Qcsk 

Vc<j  k - ]c 


RT 


(25) 


(26) 


Substitute  Eqs.  (14)  into  (26)  and  expand  in  the  spherical  co¬ 
ordinates  yields  the  modified  diffusion  equations: 
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[ar2  +F  ar  +Zjl 
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where  e\ is  the  volume  fraction  of  the  particles  type  k  in  region  j, 
and  RStkj  is  the  radius  of  the  particle  type  k.  ikj  is  the  local  transfer 
current  for  the  particles  type  k  and  is  described  by  the  Butler— 
Volmer  expression  as  follows: 


with  two  boundary  conditions  as  following 

ec^l  =  o 

dr  I  (r  =  0) 

-Mi  +z3cs,k)%^|(r=Ro4  = 


where 


(28) 


'kj 


*o  ,kj 


exp 


aa,  Li^" 
~RT~ 


-  exp 


^c,UF 

RT 


(20) 


7  =  2  Q2E 

3  “  9RT(\  -  v) 


(29) 


where  is  the  exchange  current  density  for  particles  type  k  in 
region  j  and  is  defined  as: 

=  ^Li,/<jCsurf,/<j  (cmax,/<j  —  csurf,/<j)  CLi+  (21) 


In  this  way,  to  estimate  the  stress  profile  in  the  P2D  model,  we 
only  need  to  modify  the  diffusion  equation  in  the  solid  phase  as  the 
Eq.  (27)  with  the  two  boundary  conditions  in  Eq.  (28).  When  the 
P2D  model  is  solved  as  usual  without  adding  more  variables  [39— 
41],  the  stress  profile  can  be  obtained  according  to  Eq.  (13). 


where  ku,kj  is  the  reaction  rate  constant  for  the  particles  type  k, 
Csurf ,kj  is  the  surface  concentration  of  Li+  for  the  particles  type  k, 
Cmax,kj  is  the  maximum  concentration  of  Li+  in  the  particles  type  k. 
The  over-potential  for  particle  type  k  in  region  j  is  given  by: 

Vkj  =  01  —  02  —  Ukj  ~  h<j^conc,/cj  (22) 

where  Ukj  is  open  circuit  potential  for  particle  type  k  in  region 
j ,  RCOnc,kj  is  the  contact  resistance  between  the  bulk  conductor  and 


3.  Results  and  discussion 

In  order  to  focus  on  analyzing  the  effect  of  the  adding  NCA  on 
the  LMO  material,  we  only  consider  and  discuss  the  stress  gener¬ 
ation  inside  the  LMO  particles.  First,  we  show  the  simulated  results 
for  the  pure  LMO/Carbon  cell.  And  then,  the  investigations  for  the 
cells  with  mixed  electrode  (LMO  and  NCA)  are  presented.  The 
properties  of  the  materials  and  parameters  used  in  the  following 
simulations  are  given  in  Tables  1  and  2.  Fig.  2  shows  the  simulated 


Table  1 

Properties  for  the  active  materials. 
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Parameters  Values 


LMO 

NCA 

Carbon 

Cmax  (mol  m  3) 

24161.23c 

49459.2C 

26,389d 

Ds  (m2  s-1) 

1.0e-13a 

See  Eq.  (A7)’a 

3.9e-14d 

ku  (mol-0-5  m2-5  s-1) 

5.0e-10a 

le-10a 

2.334e-llc 

R0(m) 

1.7e-6a 

2.5e-6a 

12.5e-6d 

d° 

0.19c 

0.412c 

0.65c 

Rconc  (fl  m2) 

0.05a 

0.04a 

0C 

E  (GPa) 

10b 

- 

- 

V 

0.3b 

- 

- 

Q  (mol  m-3) 

2.29e4b 

- 

- 

a  Obtained  from  Ref.  [25]. 
b  Obtained  from  Ref.  [28]. 
c  Assumed. 

d  Obtained  from  Ref.  [40]. 

results  for  a  pure  LMO/Carbon  cell  discharged  at  C/2  to  3.0  V.  The 
1  C-rate  is  given  as  10.75  A  m-2.  Fig.  2(a)  shows  the  distribution  of 
the  concentration  of  lithium  ions  on  the  surface  of  the  particles  in 
the  cathode  at  selected  times.  As  shown  in  Fig.  2(a),  the  surface 
concentration  of  particles  near  the  separator  (x  =  1 )  is  higher  than 
that  of  the  other  particles.  This  indicates  that  the  Li  ion  insertion  is 
not  homogeneous  along  the  positive  electrode  and  more  Li  ions  are 
inserted  into  the  particles  near  the  separator. 

Fig.  2(b)  shows  the  distribution  of  the  maximum  radial  stress  in 
the  particles  in  the  positive  electrode  at  the  same  selected  times  as 
in  Fig.  2(a).  The  maximum  radial  stress  is  the  largest  radial  stress 
inside  a  particle  at  the  given  time.  The  stress  profiles  presented  in 
Fig.  2(b)  show  different  trends  compared  to  the  surface  SOC  dis¬ 
tribution  presented  in  Fig.  2(a).  The  stress  in  the  particles  near  the 
separator  is  not  always  the  largest  stress.  This  is  different  from  that 
of  the  LiCo02  cathode  as  reported  by  Renganathan  et  al.  [33]. 

It  should  be  noted  that  the  maximum  radial  stress  is  located  at  the 
center  of  the  particle.  This  can  be  verified  as  follows.  Taking  the  de¬ 
rivative  of  the  radial  stress  shown  in  Eq.  (12)  with  respect  to  r  yields: 


d<rr  2QE 

( 

/  a?*) 
0 

(30) 

dr  -  3(1  -  v) 

Table  2 

Model  parameters. 

r4 

J 

Parameters 

Values 

Positive 

Separator 

Negative 

£\ 

0.41 66a 

— 

0.428b 

£ 2 

0.4a 

0.37a 

0.4b 

L(m) 

50e-6a 

50e-6a 

70e-6b 

a  (S  m-1) 

10a 

- 

100 

z 

2.89a 

3.15a 

2.5b 

aa 

0.5b 

ac 

0.5b 

t+ 

0.363a 

c°  +  (mol  m-3) 

1000a 

D  (m2  s-1) 

See  Eq.  (AS  ■ 

k  (mol  m  3) 

See  Eq.  (A4),a 

F  (C  equiv-1) 

96,487 

R  (J  mol  1  I<-1) 

8.314 

T(  K) 

298.15 

a  Obtained  from  Ref.  [25]. 
b  Assumed. 

c  Obtained  from  Ref.  [40]. 


Fig.  2.  Simulated  Li  insertion  into  LMO  particle  in  a  pure  LMO/carbon  cell  discharged 
at  C/2  to  3.0  V.  (a)  Surface  state  of  charge;  (b)  maximum  radial  stresses;  (c)  insertion 
current  density  in  the  LMO  particle  across  thickness  of  cathode  at  selected  time 
(t  =  1000,  3000,  5000,  7000,  7500  s). 
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Since  dc/dr  >  0  during  Li  insertion,  then  do>/dr  is  less  than 
zero  in  the  particle.  It  means  that  the  maximum  radial  stress  in 
a  particle  is  located  at  the  center  and  according  to  Eq.  (13)  it  is 
given  by 


max 


2QE 

9(1^7) 


Ro 

-=■  /  cr2dr  -  c 
R3nJ 


(31) 


The  two  terms  within  the  bracket  in  Eq.  (31)  are  the  average 
concentration  and  the  concentration  at  the  center  of  the  particle, 
respectively.  During  discharge  (Li  insertion  into  the  positive  elec¬ 
trode  particles),  the  concentration  at  the  center  of  the  particle  is 
the  lowest,  that  is,  the  average  concentration  is  greater  than  the 
central  concentration.  The  maximum  radial  stress  obtained  from 
Eq.  (31)  in  a  particle  is  a  positive  value,  which  means  that  the 
radial  stress  is  a  tensile  stress  during  Li  insertion.  Inversely,  during 
charge,  the  radial  stress  is  negative  and  minimum  (maximum  in 
the  absolute  value)  at  the  center  of  the  particle.  That  is,  the  radial 
stress  within  the  particle  is  a  compressive  stress  and  the 
maximum  compressive  radial  stress  is  at  the  center  of  particle.  In 
the  following  discussion,  the  maximum  radial  stress  is  the  abso¬ 
lute  value  of  the  maximum  (in  magnitude)  radial  stress  inside 
particle.  It  should  be  kept  in  mind  that  the  radial  stress  is  a  tensile 
stress  during  Li  insertion  and  is  a  comprehensive  stress  during  Li 
extraction. 

The  maximum  and  minimum  tangential  stresses  can  be  ob¬ 
tained  as  follows.  Taking  the  derivative  of  the  tangential  stress 
shown  in  Eq.  (12)with  respect  to  r  yields: 


dat 

dr 


QE 

3(1^7) 


o 


4  r4 


(32) 


Since  dc/dr  >  0  during  Li  insertion,  d<rt/dr  is  less  than  zero. 
Therefore,  the  maximum  tangential  stress  in  a  particle  is  located  at 
the  center  (r  =  0)  and  the  minimum  tangential  stress  is  located  at 
the  particle  surface  ( r  =  Ro).  The  maximum  and  the  minimum 
tangential  stresses  are  given  as  follows: 


Ro 

07, max 

 2  QE 

W0J  cr2dr-c 

~  9(1 -v) 

0 

r= 0_ 

Ro 

/TV 

QE 

W0J  cr2dr-c 

ut,min 

~  3(1— v) 

0 

r=R0. 

(33) 


By  comparing  Eqs.  (31)  and  (33),  we  can  see  that  the  maximum 
tangential  stress  is  equal  to  the  maximum  radial  stress  and  both  are 
tensile  stresses  during  the  discharging  process,  similar  to  the  nu¬ 
merical  results  reported  in  the  literature  [27].  The  minimum 
tangential  stress  is  less  than  zero  during  discharge,  because  that  the 
concentration  of  Li  ions  on  the  particle  surface  is  greater  than 
elsewhere  within  the  particle.  Since  the  minimum  tangential  stress 
is  negative,  it  is  the  maximum  compressive  stress  during  discharge. 
If  the  stress-driven  diffusion  is  ignored  in  the  diffusion  equation, 
the  analytic  solution  for  the  concentration  distribution  in  a  spher¬ 
ical  particle  can  be  obtained  as  follows:  [42]. 


cs,/<(r>0  —cs,p,  avg(t) 


hRo  ,k 

FDsk 


3_ 

10 


j'sinMlexD(iDijSi 

r  n=i  A2sin(A„)  V  R0,k  , 


(34) 


where  the  average  concentration  cS;P, avg(0  =  (3//S)  /  cr2drand 

Jo 

the  Eigenvalue  Xn  can  be  obtained  by  solving  the  Eigenequation:  tan 
(An)  =  V  The  difference  between  the  average  concentration  and  the 
central  concentration  can  be  obtained  from  Eq.  (34)  as  follows: 


Cs,p,  avg(t)  —  c|r=o  — 


0 

3  2 

FDs  k  \ 

To  2 

x  exp  | 

~Ds,k^nt 

V  Rh  . 

n  =  1  ^nSin(An) 


(35) 


The  difference  between  the  average  concentration  and  the 
surface  concentration  is  given  by: 


cs,p, avg(t)  -  c\r=R0 


kRo  ,k  (  2 

FDs,k  V10 


0°  , 

2  E  ]2eXP 

n  =  l  An 


(  ^s,/c^n  \ 

l  Rlk  )) 


(36) 


Substitution  Eqs.  (35)  and  (36)  into  (33),  we  have: 


abs(Vt?max)  —  abs((7t  min)  —  — ; 


2  QE  ikRoj< 
9(1  —v)  FDs  k 


n?i  VA"Sin(7")  '  A2 


Tlexp 


(37) 


<0 


Because  ij<  is  less  than  zero  during  discharge,  and 
En= i  (2 /AnSin(An)  +  3/^)exp(-Ds^t/JR2  k)  is  less  than  zero,  so 
Eq.  (37)  is  less  than  zero.  Eq.  (37)  shows  that  the  maximum 
compressive  tangential  stress  (o’t.min)  is  greater  than  the  maximum 
tensile  radial/tangential  stress  (o't.max/o'r.max)  during  discharge.  But, 
the  difference  shown  in  Eq.  (37)  becomes  negligible,  when  the 
concentration  gradient  is  well  developed  (t  >  (R0k/Ds,/c))-  ™S  is 
similar  to  the  numerical  results  reported  by  Christensen  et  al.  [27] 
Cheng  et  al.  have  also  derived  that  the  steady  value  of  these  two 
maximum  stresses  are  the  same  if  distribution  of  insertion  current 
density  and  stress-driven  diffusion  are  neglected. 

The  distributions  of  the  insertion  current  density  for  LMO  in  the 
positive  electrode  (z’lmo)  at  the  selected  time  were  presented  in 
Fig.  2(c).  Fig.  2(b)  and  (c)  show  that  the  distributions  of  maximum 
radial  stress  and  the  insertion  current  density  for  LMO  have  the 
same  patterns.  For  example,  at  t  =  1000  s,  both  the  radial  stress  and 
the  transfer  current  for  LMO  near  the  separator  are  maximum  in 
the  cathode.  Moreover,  it  can  be  found  that  the  maximum  radial 
stress  is  almost  proportional  to  the  insertion  current  density  at  each 
point  in  the  cathode.  The  same  phenomena  were  found  in  the  stress 
studies  with  a  single  particle  model  [26,27]  where  the  maximum 
stress  inside  a  particle  almost  linearly  increases  with  the  insertion 
current  density. 

Fig.  3  shows  the  maximum  radial  stresses  in  LMO  particles  in  a 
pure  LMO  electrode  vs.  time  at  the  interface  between  the  cathode 
and  the  separator  (Sep)  and  the  interface  between  the  current 
collector  and  the  cathode  (CC)  during  a  C/2  discharge.  The  results 
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Fig.  3.  Maximum  radial  stress  vs.  time  for  LMO  particles  at  different  positions:  current 
collector  (dash  line)  and  separator  (solid  line)  in  a  pure  LMO  electrode  discharged  at  C/ 
2  to  3.0  V. 


indicate  that  although  the  stress  in  the  particle  at  the  interface 
between  the  cathode  and  the  separator  is  not  always  the  largest 
during  the  C/2  discharge,  the  largest  insertion  stress  (around 
1000  s)  occurs  at  this  interface.  Since  the  particle  fracture  is  more 
related  to  the  largest  tensile  stress  the  particle  has  suffered, 
therefore,  the  fracture  and  loss  integrity  occur  more  likely  in  the 
particles  near  the  separator  [34].  The  fluctuation  of  the  stress 
profile  on  a  given  position  shown  in  Fig.  3  is  similar  to  Figure  1  in 
Ref.  [28  and  Figure  6  in  Ref.  [34].  As  discussed  earlier,  this  stress 
distribution  is  due  to  the  change  of  insertion  current  density  on  the 
particle  surface  over  time. 

Fig.  4  shows  the  maximum  radial  stresses  as  functions  of  applied 
current  obtained  from  two  different  models:  a  single  particle  model 
and  a  P2D  model.  The  maximum  radial  stresses  here  are  referred  to 
maximum  radial  stress  in  the  whole  electrode  during  the  discharge 
process  at  a  given  applied  current,  Japp.  In  the  single  particle  model, 
the  distribution  of  the  insertion  current  density  is  ignored  due  to 


the  assumptions  of  the  constant  electrolyte  concentration  and  the 
constant  solution  phase  potential  in  the  cell.  In  the  single  particle 
model,  we  only  solve  the  Eq.  (27)  with  two  boundary  conditions  in 
Eq.  (28).  The  insertion  current  density  i/<  on  the  particle  surface  is 
given  by: 


fappfto,LMQ 
3^1, LMO  t 


(38) 


The  stress  profile  is  calculated  based  on  the  Eq.  (13).  As  shown  in 
Fig.  4,  the  simulated  maximum  radial  stress  obtained  from  a  P2D 
model  is  higher  than  that  obtained  from  a  single  particle  model, 
especially  at  the  high  C-rates.  This  is  because  the  non-uniform 
distribution  in  the  transfer  current  in  the  cathode  is  considered  in 
a  P2D  model  as  discussed  above.  It  is  indicated  that  the  concen¬ 
tration  and  potential  gradient  in  the  solution  phase  plays  an 
important  role  in  the  mechanical  behavior  of  a  cell.  Also,  the  elec¬ 
trode  design  parameters,  such  as:  the  porosity  and  the  thickness, 
which  affect  the  transport  of  Li  ions  in  the  solution,  need  to  be 
carefully  chosen  based  on  the  application  of  the  cell.  The  effects  of 
these  parameters  on  the  cell’s  mechanical  behavior  have  been 
investigated  in  our  previous  work  [33]. 

Fig.  5  shows  the  simulated  maximum  radial  stress  for  LMO 
particles  at  selected  positions  (CC:  near  the  current  collector;  Sep: 
near  the  separator)  in  a  mixed  electrode  (LMO:NCA  =  0.67:0.33, 
volume  fraction)  which  is  discharged  at  C/2  to  3  V.  Although  the  cell 
capacity  increases  by  increasing  the  fraction  of  NCA  in  the  elec¬ 
trode,  the  1  C-rate  current  for  all  mixed  electrodes  in  this  study  is 
given  as  the  same  value  in  the  pure  LMO  cell,  which  is  10.75  A  m  2. 
The  stress  profiles  shown  in  Fig.  5  for  the  LMO  particles  in  a  mixed 
electrode  are  similar  to  those  shown  in  Fig.  3  for  a  pure  LMO 
electrode  during  a  C/2  discharge  expect  for  the  end  part  of 
discharge  (time  is  greater  than  5500  s).  As  shown  in  Fig.  5,  at  the 
end  part  of  discharge,  the  stress  in  the  LMO  particles  is  reduced 
significantly.  This  is  because  during  the  end  part  of  discharge,  most 
Li  ions  are  inserted  into  the  NCA  particles  due  to  the  OCP  difference 
between  LMO  and  NCA.  As  shown  in  Fig.  A-l  in  the  Appendix,  the 
open  circuit  potential  of  NCA  in  the  discharged  state  is  much  lower 
than  that  of  LMO. 

Fig.  6  shows  the  simulated  maximum  stresses  in  a  pure  LMO 
cathode  and  a  mixed  cathode  (LMO:NCA  =  0.67:0.33,  volume 


Fig.  4.  Maximum  radial  stress  vs.  C-rates  calculated  with  two  different  models:  porous 
electrode  model  (circle)  and  single  particle  model  (solid  line). 


Fig.  5.  Stress  profile  of  LMO  particle  at  different  position:  current  collector  (dash  line) 
and  separator  (solid  line)  in  mixed  electrode  (LMO:NCA  =  0.67:0.33,  volume  fraction) 
with  C/2  discharged  to  3.0  V. 
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Fig.  6.  Maximum  radial  stress  inside  LMO  particle  vs.  time  with  1  C  discharge/charge 
cycling  of  mixed  electrode  (black,  LMO:NCA  =  0.67:0.33,  volume  fraction)  or  pure 
electrode  (green)/carbon  full  cell  between  3.0  V  and  4.3  V.  The  ellipse  indicated  the 
stress  reduction  region  in  the  mixed  electrode  compared  to  pure  electrode  during 
cycling.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  this  article.) 


fraction)  as  functions  of  time  during  a  1  C  discharge/charge  cycle. 
The  black  line  represents  the  stress  profile  obtained  from  a  blended 
material  cathode;  the  green  (in  web  version)  line  is  maximum 
stress  profile  for  the  pure  LMO  cathode.  As  shown  in  Fig.  6  (red 
ellipse  (in  web  version)),  the  stress  of  LMO  particle  in  the  dis¬ 
charged  state  is  reduced  due  to  mixed  with  the  NCA  particles 
compared  that  obtained  from  a  cell  with  a  pure  LMO  cathode. 

Some  experimental  work  has  shown  that  the  dissolution  of  Mn 
from  spinel  LMO  is  a  serious  process  in  a  discharged  state  during 
cell  cycling  [43-45].  Mn  dissolution  into  the  electrolyte  due  to  acid 
attack  and  a  disproportionation  reaction  at  the  particle  surface  is 
most  likely  to  happen  with  more  Mn3+  inside  particle  when  the  cell 
is  in  the  discharged  state.  If  there  is  a  protective  layer  on  the  particle 
surface  such  as  a  coating  layer  or  a  SEI  layer,  the  Mn  dissolution  can 


be  reduced  [9  .  As  determined  by  experiments,  the  inactive  layer  is 
always  formed  on  the  LMO  particle  surface  due  to  the  Mn  disso¬ 
lution  and  the  electrolyte  oxidation  [46,47].  However,  unlike  the 
bulk  LMO  particles  which  have  been  treated  at  high  temperature, 
the  newly  formed  layer  may  be  weak  and  be  easily  destroyed  due  to 
the  stress  as  shown  in  Fig.  7.  Therefore,  the  decrease  in  the  stress  in 
the  discharged  state  for  LMO  particles  in  the  mixed  cathode  may 
suppress  the  fracture  of  the  SEI  layer.  Therefore,  the  existence  of  the 
layer  in  the  discharge  state  can  effectively  protect  the  LMO  surface 
from  the  acid  attack  and  reduce  the  Mn  dissolution.  Moreover,  if  the 
stress  on  the  particle  surface  contributed  to  the  Mn  dissolution 
directly  due  to  the  increase  of  strain  energy,  the  Mn  dissolution 
would  be  suppressed  when  the  stress  is  reduced.  As  discussed 
above,  the  reduction  of  the  stress  in  a  discharged  state  can  explain 
the  experimental  observations  that  the  performance  of  LMO  based 
electrode  are  improved  by  adding  NCA  [16]. 

The  maximum  radial  stress  of  the  cell  with  mixed  electrode 
(LMO:NCA  =  0.67:0.33,  volume  fraction)  during  a  constant  current 
discharge/charge  cycles  (CC-CC  protocol)  at  1  C  is  shown  in  Fig.  8. 
The  results  indicate  that  the  stresses  in  the  latter  cycles  are  similar 
to  that  in  the  first  cycle  in  this  study.  This  is  because  the  current  is 
low,  the  effect  of  previous  Li  insertion  and  extraction  on  the  stress 
vanishes  within  a  short  timescale  [27].  The  maximum  radial  stress 
in  a  cell  cycled  with  the  CC  discharge-CC  charge  protocol  is  pre¬ 
sented  in  Fig.  9.  This  protocol  consists  of  a  constant  current 
discharge  at  1  C  to  3.0  V,  followed  by  a  constant  current  charge  at 
1  C  to  a  selected  cutoff  charge  voltage,  and  then  a  constant  voltage 
charge  until  the  current  drops  to  the  cutoff  current.  Similarly  to  the 
CC-CC  protocol  cycle,  there  is  no  stress  accumulation  during  the 
cell  cycling  in  this  study.  Also,  the  simulation  shows  that  the  stress 
generation  during  the  CV  charge  with  a  cutoff  voltage  of  4.3  V  is 
similar  to  that  with  a  cutoff  voltage  of  4.2  V.  The  simulated 
maximum  stress  profiles  for  the  mixed  LMO  and  NCA  electrodes 
with  four  different  volume  ratios  are  shown  in  Fig.  10.  The  following 
volume  ratios  between  LMO  and  NCA  are  selected  in  this  study: 
LMO: NCA  =  0.9:0.1,  0.67:0.33,  0.5:0.5,  0.33:  0.67.  The  cells  are  dis¬ 
charged  at  3  C  to  3.0  V.  As  expected,  the  reduced  stress  region  in¬ 
creases  with  increasing  the  volume  fraction  of  NCA.  But  it  is  also 
observed  from  Fig.  10  that  the  maximum  stress  during  discharge  (as 
indicated  by  the  dash  arrow)  increases  by  adding  more  NCA  into  the 
mixed  electrode.  Therefore,  the  tradeoff  between  the  decrease  in 
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Fig.  7.  Schematic  of  failure  that  may  occur  in  LMO  particle  surface  due  to  Li  insertion/extraction:  (a)  Li  insertion  (discharge):  (b)  Li  extraction  (charge). 


Y.  Dai  et  al.  /  Journal  of  Power  Sources  247  (2014)  365-376 


373 


Fig.  8.  Simulated  stress  generation  inside  LMO  particle  with  1  C  discharge/charge 
cycling  of  mixed  electrode  (LMO:NCA  =  0.67:0.33,  volume  fraction)  between  3.0  V  and 
4.3  V.  The  dash  line  is  the  applied  current,  corresponding  to  right  axis.  Positive  current 
is  charged  process,  and  negative  current  is  discharged  process. 


the  stress  in  the  end  part  of  discharge  and  the  increase  in  the 
maximum  stress  has  to  be  considered  when  determining  the  vol¬ 
ume  fraction  of  NCA  in  the  mixed  electrode. 

4.  Conclusions 

A  mathematical  model  is  developed  in  this  study  to  understand 
the  stress  generation  inside  the  LMO  particle  during  discharge/ 
charge  in  a  Li-ion  battery  with  a  pure  LMO  or  a  mixed  electrode 
(LMO  and  NCA)  cathode.  Simulation  results  for  the  pure  LMO/car- 
bon  full  cell  show  that  the  stress  generation  inside  LMO  particle  is 
non-uniform  within  the  positive  electrode.  The  results  also  indicate 
that  the  fracture  or  loss  of  particle  integrity  potentially  occurs  in  the 
particles  near  the  separator.  The  stress  profiles  obtained  from  a 
single  particle  model  and  a  porous  electrode  model  show  a  sig¬ 
nificant  difference  for  the  high  C  rate  discharge.  This  indicates  that 
the  single  particle  model  cannot  be  used  because  of  the  non- 


Fig.  9.  Simulated  stress  generation  inside  LMO  particle  in  mixed  electrode 
(LMO:NCA  =  0.67:0.33,  volume  fraction)  with  1  C  discharge/1  C  charge  to  cut-off 
voltage  and  constant  current  charge  to  cut-off  current.  Two  different  cut-off  charge 
voltages,  4.2  V  (solid  line)  and  4.3  V  (dash  line)  are  applied.  The  arrow  in  the  picture 
indicates  a  new  cycle. 


Fig.  10.  Simulation  stress  vs.  time,  of  different  blend  ratio  mixed  electrode  with  3  C 
discharged  to  3.0  V.  0.33  -  red,  0.5  -  black,  0.67  -  purple,  0.9  -  green  of  LMO  volume 
fraction  in  mixed  electrode.  The  dash  arrow  indicates  the  largest  stress  during  the 
whole  discharge  process.  (For  interpretation  of  the  references  to  colour  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


uniform  distributions  in  the  concentration  and  the  potential  in 
the  electrolyte  in  the  high  current  rate  studies.  Simulation  results 
for  a  mixed  electrode  show  that  the  stress  generated  in  the  LMO 
particles  is  significantly  reduced  at  the  end  of  discharge  in  a  mixed 
LMO  and  NCA  electrode  due  to  the  OCP  difference  between  the  LMO 
and  NCA.  This  stress  reduction  might  partially  explain  why  exper¬ 
iment  investigations  show  that  the  performance  of  the  LMO  based 
electrodes  can  be  improved  by  adding  NCA.  Simulations  for  mixed 
electrodes  with  different  blend  ratios  show  that  the  reduced  stress 
region  increases  by  increasing  the  NCA  ratio. 

Appendix 

The  effective  ionic  conductivity  and  diffusion  coefficient  in  the 
binary  electrolyte  are  determined  by  the  following  equations, 
respectively:  [36]. 

Keff J  =  j  =  p,s,n  (Al) 

Tj 

Deffj  =  —A  j  =  P,s,n  (A2) 

Tj 

where  z j  is  the  tortuosity  of  the  porous  electrode. 

The  concentration  dependent  ionic  conductivity  and  the  diffu¬ 
sion  coefficient  in  the  binary  electrolyte  are  given  by:  [36]. 


(i  13/1  CLi+  \ 

- ,  '  .  1000  . - -4  +  O.I  (A3) 

1  +  0-2  (nftio)  +  0-08  (tsqo)  J 

D  =  6.5  X  1010exp(  -  0.7^)  (A4) 

The  effective  conductivity  in  the  solid  phase  is  defined  as 
ffeff j  =  <Tj,  j  =  p,n  (A5) 
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The  ionic  activity  coefficient,/,  is  given  as  following:  [36]. 


i+/3L=i 

9lncLi+ 


1000 


-1.0178 


2  (1  +0.983 1(^)  '  ) 
0.9831 


(too)  1  +  0.9831  (^Sjj 


0.5 


1.5842 


(A6) 


The  Li  diffusion  coefficient  in  NCA  particle  [36]  is  given  as 
follows: 


Ds.nca  =  3  x  10  15((1  +  tan  h(  -  2O(0NCA  -  0.73)))  +  0.02) 

(A7) 

The  open  circuit  potentials  for  the  LiMn2C>4  cathode  as  functions 
of  state  of  charge  are  given  by:  48]. 


UNCA  =  8.535  -  17.O590NCA  +  21.0380^  -  9.1530^ 

+  9.875(0nca  -  0.700)3  -  537O.872(0NCA  -  0.935)3; 

(A13) 


0.959  <  ^nca  —  0.980, 

UNCA  =  8.535  -17.O590NCA +21 .0380^* -9.1530]^ 

+  9.875(0nca  -  0.700)3  -  537O.872(0nca  -  0.935)3 
-  47690. 3O4(0nca  -  0.959)3. 

(A14) 

And  the  open  circuit  potential  for  the  carbon  anode  (fit  to 
experimental  data)  is  as  follows: 

0.001  <  0C  <  0.0109, 


ULM0  =4.199  +  0.05661  tan  fi(-14.5550LMO  +  8.609) 

;(-  1 


-0.0275 


\  0.492 


1.901 


\(0.998  -  f7LMOj 

-  O.157exp(-O.O4740LMO) 

-  0.810exp(  -  4O(0LMo  -  0.134)). 


(A8) 


The  open  circuit  potential  for  NCA  (fit  to  experimental  data)  is 
given  by: 


Uc  =  0.113  -  0.0208  tan  h  (15.O640C  -  8.199) 

-  2.435(0C  -  0.440)3  +  65.394(0C  -  0.154)3 

-  96O.3O7(0C  -  0.0897)3  -  1.006  x  1O7(0C  -  0.0109)3; 

(A15) 

0.0109  <  0C  <  0.0897, 


0.36  <  0NCA  <  0.410, 

UNCA  =8.535-  17.0590^  +  21.0380^-9  1530^ 

+  9.875(0nca  -  0.700)3  -  2.176(0nca  -  0.550)3 
—  1331 .866(0NCA  —  0.41 0)3; 

(A9) 

0.410  <  0NCA  <  0.55, 

UNCA  =  8.535  -  17.O590NCA  +  21.0380^  -  9.1530nca 
+  9.875(0nca  -  0.700)3  -  2.176(0nca  -  0.550)3; 

(A10) 


Uc  =  0.113  -  0.0208  tan  h  (15.O640C  -  8.199) 

-  2.435(0C  -  0.440)3  +  65.394(0C  -  0.154)3 

-  96O.3O7(0C  -  0.0897)3;  (A16) 

0.0897  <  0C  <  0.154, 

Uc  =  0.113  -  0.0208tan  h  (15.O640C  -  8.199) 

-  2.435(0C  -  0.440)3  +  65.394(0C  -  0.154)3;  (A17) 

0.0897  <  0C  <  0.154, 

Uc  =  0.113 -0.0208tanh  (15.O640C -8.199) 

-  2.435(0C  -  0.440)3  +  65.394(0C  -  0.154)3;  (A18) 


0.55  '  ■  0jv|ca  '  0.7, 


0.154  <  0C  <  0.440, 


Unca  =8.535  -  17.0590^  +  21.0380^  -  9.1530^ 
+  9.875(0nca  —  0.700)3; 

0.7  <  ^nca  —  0.935, 

UNCA  =  8.535  -  17.O590NCA  +  21.0380^  -  9.1530^ 


Uc  =  0.113  -  0.0208tan  h  (15.O640C  -  8.199) 

-  2.435(0C  -  0.440)3; 

(AH)  0.440  <  0C  <  0.854, 

Uc  =  0.113-  0.0208tan  h  ( 1 5 ,0640C  -8.1 99) ; 
0.854  <  0C  <  0.92, 

(A12) 

Uc  =  0.113  -  0.0208tan  h  (15.064 0C  -  8.199) 

-  252.7O7(0C  -  0.854)3; 


(A19) 


(A20) 


0.935  ■  0nca  —  0.959, 


(A21) 
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where  the  SOC,  6,  is  defined  by: 

0k  =  5?H k  =  LMO,  NCA,  C  (A22) 

^max,/c 


SOC 

Fig.  A-l.  The  open  circuit  potential  profiles  for  (a)  LMO;  (b)  NCA. 


List  of  symbols 

a/<  specific  surface  area  of  particle  type  k,  m2  m-3 

bruggj  Bruggeman  coefficient  in  region  j  (j  =  n,  s,  p) 
cLi+  concentration  of  Li+  in  the  solution  phase,  mol  m-3 

cs,k  concentration  of  lithium  ions  in  the  particle  type  k, 

mol  m  3 

c°s,k  initial  concentration  of  lithium  ions  in  the  particle  type  k, 
mol  m  3 

cs, max, k  maximum  concentration  of  lithium  ions  in  the  particle 
type  k,  mol  m-3 

A cs  concentration  change  from  initial  value,  mol  m-3 

D  Li+  diffusion  coefficient  in  the  electrolyte,  m2  s_1 

Ds<k  lithium  ion  diffusion  coefficient  in  the  particle  type  k, 
m2  s_1 

E  Young’s  modulus,  GPa 

/  ionic  activity  coefficient 

F  Faraday’s  constant,  96,487  C  equiv-1 

/app  applied  current  density,  A  m-2 


io,k  exchange  current  density,  A  m-2 

ik  transfer  current  on  the  surface  of  particle  type  k,  A  m”2 

Ij  total  current  density,  A  m~3 

1<u,k  Li  intercalation/deintercalation  reaction  rate  constant  for 
particle  type  k,  A  m25  mol- 15 
Lp  thickness  of  positive  electrode,  m 

Ls  thickness  of  separator,  m 

Ln  thickness  of  negative  electrode,  m 

n  negative  electrode 

p  positive  electrode 

r  radial  coordinate,  m 

R  universal  gas  constant,  J  mol-1  K-1 

Ro,k  radius  of  the  particle  type  k,  m 

Fconc  contact  resistance  on  the  surface  of  particle  type  /<,  Q  m2 
s  separator 

t  time,  s 

t+  Li+  transference  number  in  the  electrolyte 

T  environment  temperature,  K 

AT  temperature  change  from  initial  value,  K 

u  displacement,  m 

Uk  equilibrium  potential  of  Li  interaction/deinteraction 
reaction  for  particle  type  /<,  V 
x  spatial  coordinate,  m 

Zi,  Z2  integration  constant 


Greek  letters 

a  thermal  expansion  coefficient 

aa  anodic  transfer  coefficient 

ac  cathodic  transfer  coefficient 

e\ tk  volume  fraction  of  the  particle  type  k  in  the  electrode 

€2j  porosity  in  region  j  ( j  =  n,  s,  p) 

eij  strain 

6  dimensionless  concentration  of  lithium  ions  in  the 

particle  {6  =  cSiSUrf/cs,max) 

6 0  initial  dimensionless  concentration  of  lithium  ions  in  the 

particle  (6  —  cStSurf/ cSimax) 

k  ionic  conductivity  of  the  electrolyte,  S  m-1 

Ke ff  effective  ionic  conductivity  of  the  electrolyte,  S  m_1 

a  electronic  conductivity  of  the  solid  phase,  S  m-1 

(7eff  effective  electronic  conductivity  of  the  solid  phase,  S  itT1 

( Tij  stress,  Pa 

solid  phase  potential,  V 
02  solution  phase  potential,  V 

r]k  over-potential,  V 

£  dummy  integration  variable 

v  Poisson’s  ratio 

6y  Kronecker  delta 

t  tortuosity 

Q  partial  molar  volume,  m3  mol-1 


Subscript 

h  hydrostatic 

r  radial  direction 

t  tangential  direction 

LMO  LiMn204 

NCA  LiNio.8Coo.15Alo.05O2 

C  carbon 
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